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Abstract 

We investigate the Yang-Baxter algebra for U(l) invariant three-state vertex models whose 
Boltzmann weights configurations break explicitly the parity-time reversal symmetry. We 
uncover two families of regular Lax operators with nineteen non-null weights which ultimately 
sit on algebraic plane curves with genus five. We argue that these curves admit degree two 
morphisms onto elliptic curves and thus they are bielliptic. The associated R-matrices are non- 
additive in the spectral parameters and it has been checked that they satisfy the Yang-Baxter 
equation. The respective integrable quantum spin-1 Hamiltonians are exhibited. 
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1 Introduction 



Over the past decades we have witnessed the importance played by vertex models in the de- 
velopment of the theory of integrable systems in two spatial dimensions [TJ. Let us recall some of 
the basic notions about this classical lattice model of statistical mechanics. For simplicity consider 
a square lattice of size N x N with periodic boundary conditions on both horizontal and vertical 
directions. The statistical configurations are specified by assigning to each set of four edges meeting 
at a given lattice site the spins variables a, /3, 7 and 5. Here we assume that these state variables 
take values on a finite subset V of the integers, that is V = {1, • • • ,q}. To a given vertex of 
the lattice we assign a Boltzmann weight W^'t representing the energy of the corresponding spins 
configurations. In figure [I] we illustrated our notation for the weight indices. 
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Figure 1: The Boltzmann weight of the vertex model on a square lattice. 

One feature of the vertex models is that its basic properties can be formulated with the help 
of a beneath tensor structure formally represented by the product Vo <S> 0. The horizontal 
degrees of freedom are encoded in the term Vo which is often called auxiliary space. The second 
term Vk stands for the vertical degrees of freedom associated to each k-th site of a one dimensional 
lattice of size N playing the role of the quantum space of a g-state spin chain. It turns out that 
the Boltzmann weights can be organized in terms of a local matrix acting on Vq <S> denominated 
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Lax operator, 

L fc (w) = £ W# egjseg, for fc = 1, • • ■ , N, (1) 

a, /3, 7, <5=1 

where denotes g x g Weyl matrices acting on the space Vj with j — 0, • • • , N. From now on we 
shall refer to the set of nonzero weights Wjt using the symbol "w" . 

The respective row-to- row transfer matrix T(w) can then be written as the trace over the 
auxiliary space of an ordered product of Lax operators, namely 

T(w) = Tr Vo [L N (w)L N _i(w) • • • L 2 (w)Li(w)] . (2) 

A relevant family of vertex models are those whose Lax operators are invariant by a single U(l) 
symmetry. This invariance implies that many of the weights are zero depending on whether or not 
the respective indices satisfy the so-called ice condition, 

•W^ = 0, for « + /3^ 7 + 5 

•WJJ^O, for a + /3 = j + 5. (3) 

Up to the present the known realizations of vertex models satisfying the above rule have the 
corresponding weights parameterized in terms of trigonometric functions. The typical examples are 
the vertex models associated to the solutions of the Yang-Baxter equation based on the quantum 
U[SU(2)]§ algebra either for generic values of the deformation parameter (HIE] or when it takes 
values on the roots of unity OE]. The current results in the literature suggest that to obtain 
integrable vertex models with weights not uniformized by rational functions one has to consider 
statistical configurations that violate the U(l) symmetry. For example, these are the cases of 
certain generalizations of the eight- vertex model [7] and the celebrated chiral Potts model [HIE] 
having both an underlying Z g symmetry. 

On the other hand, it has been shown that the transfer matrix of U(l) invariant vertex models 
can be diagonalizable by the algebraic Bethe ansatz for arbitrary Lax operators without reference 
to a given specific parameterization of the weights pi)] . In this work it was not needed to make any 
assumption on the dependence of the spectral parameters entering the corresponding R-matrix to 
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build up the transfer matrix eigenvectors. In any way this algebraic approach forbids the existence 
of vertex models having both the U(l) invariance and Boltzmann weights sitting on algebraic 
varieties which can not be rationally uniformized. In fact, the generality of a number of weights 
identities derived in [10] make it hard to believe that they are only realized in terms of trigonometric 
functions. 

Of course, for such class of models, irrationality of weights can only emerge when the number 
of states q is larger than two. This is because the model with q = 2 corresponds to the asymmetric 
six- vertex model whose weights are known to be rationally parameterized, see for example [II] . For 
q > 2, however, the structure of the functional relations derived from the Yang-Baxter algebra is 
very different from that satisfied by weights of an arbitrary six-vertex model. Would they be so 
stringent to always drive us to exactly solvable vertex models with trigonometric weights?. Here 
we investigate this question in the most simple case where non-rational weights can not be rule 
out: the three-state U(l) vertex model. Another relevant motivation to study these kind of systems 
comes from the existence of concrete exactly solvable spin-1 quantum chains discovered by Alcaraz 
and Bariev within the coordinate Bethe ansatz method [12]. The fact that their Hamiltonian for 
general couplings can not be derived in terms of an additive R-matrix suggests that non rational 
three-state vertex models should indeed exist. 

Our study of Yang-Baxter algebra for U(l) three-state vertex models will lead us to develop a 
strategy to deal with a problem involving a large number of functional relations constraining the 
Boltzmann weights. The basic guidelines of our approach is somehow general and can in principle 
be used to study more complicated vertex models. It turns out that we are able to uncover two 
families of integrable U(l) three-state vertex models with weights lying on non-rational manifolds. 
In fact, we shall argue that their Boltzmann weights are ultimately constrained by algebraic plane 
curves of genus five. We remark that the respective quantum spin-1 chains extend in a substantial 
way the previous integrable Hamiltonian found in reference [12] . 

We have organized this work as follows. In next section we describe some basic properties of 
the polynomial relations coming from the Yang-Baxter algebra. We have endeavored to make it 
self contained when useful mathematical notions of algebraic geometry are used. In section [3] we 
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describe the main structure of the U(l) invariant three-state vertex model to be studied in this 
paper. The solution of the respective functional relations is detailed in Section 0] and we have been 
able to uncover two families of integrable three-state vertex models. In Section [5] we investigate 
the geometrical properties underlying the integrability of these vertex models and for one of the 
families this forces us to analyze the problem of the intersection of two projective surfaces. We 
found that the properties of the underlying algebraic manifolds are related to that of genus five 
bielliptic curves. We have summarized the main results for the Lax operators and the respective 
R-matrix in Section [6j We compute the expressions of the corresponding exactly solvable spin-1 
chains and show that they contain as particular case the spin-1 Hamiltonian found previously by 
Alcaraz and Bariev [T2]. Our conclusions are presented in Section [7] and in four Appendices we 
describe a number of technical details complementing the discussions of the main text. 



2 Integrability Conditions 

In general, a lattice model of statistical mechanics in two-dimensions is considered integrable 
when the corresponding transfer matrix can be embedded into a family of pairwise commuting 
operators [I], 

T(w'),T(w") 

where w' and w" represent two different sets of weights. 

A sufficient condition for commuting transfer matrices was originally introduced by Baxter in 
his analysis of the eight- vertex model [13]. This condition requires the existence of an invertible 
R-matrix which together with the Lax operators should satisfy the Yang-Baxter algebra, 

R(w',w")[L fc (w')®I 3 ][I g ®L fc (w")] = [I,<8>L fc (w")][L fe (w') ®LjR(w',w"), (5) 

where l q denotes the q x q identity and R(w', w") is a q 2 x q 2 matrix acting on the tensor product 
V <8> V . 

The commutation relation of two distinct transfer matrices (jl]) is obviously not affected when 
their weights are multiplied by two independent nonzero scalar factors. This means that the 
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functional equations coming from Yang-Baxter algebra are expected to be homogeneous separately 
in each of the sets w' and w" of weights. To be more precise writing Fj(w', w") to denote a given 
polynomial derived from Eq.((5]) we then, in general, can state that, 

F j (X l w',X 2 w") = Af^F^wV), V Ai, A 2 ^ 0, (6) 

where D\ and D 2 define the bidegree of the bihomogeneous polynomial Fj(w',w"). 

In this paper we will also assume that the R-matrix satisfies the standard unitarity condition, 

R(w', w")PR(w", w')V = p(w', w")I q <g) I q (7) 

where V denotes the permutator operator acting on a g 2 -dimensional space and p(w', w") represents 
an overall normalization. 

The above assumption is motived by the fact that unitarity property ([7]) assures us from the very 
beginning that the R-matrix has an inverse. Recall that unitarity has also been relevant in providing 
us a number of identities that were essential for the algebraic diagonalization of the transfer matrix 
of the U(l) invariant vertex models [in]. In addition, we shall show that the unitarity property of 
the R-matrix imposes an important restriction on the structure of the polynomials Fj(w',w"). In 
order to see that we multiply the left and right sides of Eq. (J5]) by the inverse of the R-matrix and 
with the help of Eq.([7]) we obtain, 

[Lfc(w') ®y[I ? «)L fc (w")]PR(w",w')P = PR(w",w')P[I,®L fc (w")][L fc (w')®IJ. (8) 

We now apply the permutator on both sides of Eq.© as well as we insert the identity V 2 = l q <S)l g 
in the middle of the brackets to permute the Lax operators. As a result we can derive the following 
relation, 

[Iq®L fe (w)][L fc (w )<g>LjR(w ,w) = R(w , w)[Lfc(w ) <g) I ? ][I 9 <8> L fc (w )]. (9) 

Inspecting Eqs. (!5|9|) we see that their left and right sides are related once we interchange the 
weights, that is w' -h- w". This means that the polynomial equations coming from the Yang-Baxter 
algebra are expected to be anti-symmetrical upon the exchange of weights label, that is, 

F^w', w") + F^w", w') = 0. (10) 
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We stress that such simple consequence of the unitarity of the R-matrix is going to play an 
important role to help us disentangle involved high degree functional relations on the Boltzmann 
weights. It is however fortunate that in many instances of our analysis this fact will come out 
naturally since we will be able to write the polynomials in the following particular anti-symmetrical 
form, 

F,(w', w") = H,(w')G,(w") - H.VXV), (11) 

where Hj(w) and Gj(w) are irreducible homogeneous polynomials with the same degree D. We 
then say that Fj(w',w") is an irreducible bihomogeneous polynomial with bidegree (D,D). 

We see that polynomials of the form (fTTj) vanish trivially when we consider the limit w' — > w" 
which is a desirable property since certainly the transfer matrix commutes with itself. We next note 
that such bihomogeneous polynomials always admit an special solution in which the distinct group 
of weights w' and w" are decoupled from each other. This solution bears some resemblance with 
the method of separation of variables used to solve the classical dynamics by the Hamilton- Jacobi 
theory and partial differential equations of mathematical physics. It can be written as follows, 

where the parameter Aj is considered a free constant. 

It turns out that such particular solution to Eq.f fTT]) has a very clear meaning in the realm of alge- 
braic geometry. This discussion permits us to introduce the appropriate mathematical terminology 
making presentation self contained. We start by recalling that the zero locus of the bihomogeneous 
polynomial Fj(w',w") is known to produce a well defined algebraic variety [13] . In order to 
define this mathematical object let us denote the set of weights by the elements Uq, • • • , w m repre- 
senting the coordinates of a projective space CP m over the complex field. The algebraic variety Xj 
is a closed subset of the product of such two projective spaces which formally can be represented 
as, 

X, = {[u>' :...: u'j x ^ : . . . : J' m \ e CP m x CP m |H,(4, . . . , ujO^, . . . , O 

- R j (u)o,...,u)'^)G j (u)' ,...,u m ) = 0}, (13) 
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where [ojq : . . . : oj m ] denotes a point in the projective space CP m by which we mean the line 
spanned by the vector (oj , . . . , u) m ) G C m+1 where the origin is omitted. 

By the same token, the polynomials originated from the special solution (112j) can also be used 
to define an underlying subvariety Yj C Xj. This subvariety is in fact described by the product of 
two identical algebraic sets since the corresponding polynomials do not mix distinct weights labels. 
We then are able to write Yj = Y(Aj-) x Y(Aj) where the component Y(Aj) is defined by, 

Y(A i ) = {[tu : . . • : u m ] G CP m |Hj(u;o, ■ ■ ■ , u m ) - AjG^o, ...,co m ) = 0} (14) 

We now can show that the particular solution ( 1T2|) gives rise to a divisor on the original variety 
Xj. To this end we recall that one basic invariant of any variety is its dimension which here can 
be determined using the standard result that an irreducible hypersurface S(wq, . . . ,w m ) G CP m 
has dimension dimS = m — 1 [H]. From this result it follows that the variety Xj has dimension 
dimXj = 2m — 1 while the dimension of the subvariety Yj is dirnY^ = 2m — 2 since they are 
generated by irreducible polynomials. From the fact that dimX^ — dimY^ = 1 and observing that 
the intersection multiplicity of Yj at Xj is also 1 we then conclude that Yj is in fact a prime divisor 
element on Xj. This means that by varying the parameter Aj we are able to foliate the variety Xj 
through submanifolds of codimension 1 whose fibers are determined by the variety Y(Aj). 

We shall see that the prime divisors associated to polynomials with the structure (TTTT) are 
precisely the fundamental building blocks of a vertex model with commuting transfer matrix. It 
is the intersection of a collection of such divisors that ultimately is going to dictate the algebraic 
variety in which the Boltzmann weights are lying on. This procedure assures us the existence of two 
independent transfer matrices that are sited on the same algebraic manifold and thus of a single 
family of Lax operators. 

3 Three-state Vertex Model 

We now turn our attention to the presentation of the specific U(l) three-state vertex which 
we intend to investigated in this paper. From the ice-rule ([3]) it follows that this type of model 
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can have at most nineteen different Boltzmann weights. This space of parameters can be reduced 
once we consider typical symmetries of the weights when they are viewed as (1 + l)-dimensional 
scattering amplitudes pj)]. These invariances are defined as follows, 



►Parity Reversal : W 



7,5 



w 



5,7 
/3,a> 



•Time Reversal : W 



7,5 



•Charge Conjugation : W 



W 

7,5 



a,/3 
7,5 i 



(15) 
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3+1 — f,q+l—6 



According to the recent work [16] nineteen vertex models invariant by the combined action 
of parity and time reversal symmetries have always rational weights. This means that we have 
to consider vertex models whose statistical configurations do not preserve the PT transformation. 
From Eqs.(fl5]) we see that this is achieved when the respective weights fulfill one of the following 
inequalities, 
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(16) 



One way to assure the breaking of the PT symmetry is by means the diagonal weights 
since the off-diagonal ones can in principle be modified with the help of gauge transformations. 
In this case, broken PT invariance is not completely incompatible with the preservation of charge 
conjugation which in turn permits us to work with a smaller number of distinct weights. Considering 
that charge symmetry is preserved at least by the diagonal weights our starting ansatz for the Lax 
operator is, 
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(17) 



where a, b, b, c, d, J, /, g, h and denote ten distinct weights of the set w. We see that the PT 
invariance is only broken by way of the diagonal weights b and b. 

We now consider an arbitrary R-matrix and substitute it together with the above ansatz for the 
Lax operators in the Yang-Baxter algebra (jSJ). It is not difficult to see that the underlying U(l) 
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invariance of the Lax operators imposes us severe constraints on the R-matrix. Under the mild 
assumption that some of the weights of the Lax operators are not trivially related we find that 
R-matrix elements have also to satisfy the ice-rule ([3]). This motivates us to choose the R-matrix 
with the same structure of the Lax operators, namely 



R(w , w 
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(18) 



where bold letters are used to distinguish the R-matrix elements from the Boltzmann weights. 

At this point we have the basic ingredients to study the possible solutions of the Yang-Baxter 
algebra ([5]). We shall tackle this problem using the following systematic strategy. We start by elim- 
inating the elements of the R-matrix since the main purpose is the determination of weights fixing 
the Lax operators. To this end we search for suitable functional relations that built a consistent 
linear system of homogeneous equations for a particular chosen subset of R-matrix entries. The 
vanishing of respective determinant is going to lead us in most cases to polynomials having the 
anti-symmetrical structure ffTTj) . This makes it possible to define the associated divisors ( lT2j) and 
as a result the freedom of a number of free parameters Aj. In the situation of functional relations 
that can not be written directly in the special form ({TTj) we impose that their polynomials should 
satisfy the anti-symmetric property (TTU]) . This idea is crucial to solve very involved functional 
equations resulting from many nested steps. It turns out that we always will be able to implement 
this property at the expense of imposing constraints among the free parameters Aj. As a result 
either the corresponding polynomial vanishes directly or it can be brought into the suitable form 
f fTTj) . After all that, we still have to perform the intersection of some basic divisors which is going 
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to lead us to the main algebraic manifold for the Boltzmann weights. As a byproduct we are able 
to determine the matrix elements of the R-matrix in terms of few Lax operators weights. In next 
sections we show how to carry out all these steps in practice. 

4 The Functional Relations 

The functional equations constraining the entries of the R-matrix and the Boltzmann weights are 
derived by substituting our proposals ( ITTffTS]) in the Yang-Baxter algebra. We find that there exists 
fifty-four independent equations which are best subdivided in terms of their number of distinct 
terms. In Table [1] we summarize this classification which ranges from relations having only two 
terms to those with the maximum number of five elements. 



Number of Equations 


Number of Terms 


2 


two 


15 


three 


25 


four 


12 


five 



Table 1: The number of distinct functional relations versus their respective number of terms. 



In general, the technical difficulties in dealing with the solution of the functional equations 
increase with their number of terms due to the presence of many free variables. However, we shall 
see that the analysis of the simplest relations with two terms will result in a reasonable decrease 
of the number of independent equations having three elements. This simplification is important to 
make possible the solution of such functional equations by the elimination method. In addition to 
that we will need to analyze only a small number of equations with four terms, as compared to the 
available set of Table [H to decide about the integrability of the vertex model we have started with. 
After this study we are left to verify that all the remaining polynomial equations coming from the 
Yang-Baxter relation are satisfied. It turns out that this task can be performed algebraically with 
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the help of computer algebra system. 



4.1 Two Terms Equations 

The functional equations having two terms are given as follows, 



c(d' d" — d'ct) = 0, 
(dd! - dtT)c" = 0. (19) 

We see that the above equations are already in the convenient form ffTTl) since we are disregarding 
possible solutions with zero weights and R-matrix amplitudes. Clearly, the corresponding divisors 
fix the following ratios among the variables, 

d d , . 

d = d = A °- < 20 > 

Here we are tacitly assuming that A ^ since otherwise we would have to set the weight d to 
zero. Our main interest is to consider genuine nineteen vertex models and thus all the Lax operator 
weights must be non-null. We now turn to investigate more complicated functional relations having 
three and four terms. 
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4.2 Three Terms Equations 

By considering the ratios ( 12 Up the number of independent relations with three elements decrease 
from the original fifteen to only nine functional equations. Their expressions are given by, 



/ // — / // / // 

ac a — be b — ca c = 0, 


(21) 


in in i a 
ccb + ha c -aoc =0, 


(22) 


ad'b" — cb' d" — be?'/" = 0, 


(23) 


, i /I r i ji j (.a 
bo a — a/ a + cd j = 0, 


(24) 


hd'a" — fd'b" — db'c" = 0, 


(25) 


ddb" + fb'c" - bf'c" = 0, 


(26) 


— / // /— // — / // 

cb a —cab — be c =0, 


(27) 


cf %" _ bdV - cb'f" = 0, 


(28) 


df'a -db'b" -fd'c" =0. 


(29) 



We first observe that Eqs. (12 11129]) are not invariant when we exchange the R- matrix elements 
with the corresponding double primed Lax operator weights. The main reason for the absence of 
this invariance is because we are assuming broken PT symmetry that is b ^ b and b ^ b. Therefore, 
there exists a concrete possibility that the R-matrix and the Lax operators may be sited in two 
distinct algebraic varieties that are not isomorphic. This is already an indication we have some 
chance to obtain an integrable vertex model whose weights are not trigonometric. 

We now consider the solution of Eqs.( l211 l 2"9~|) as a system of homogeneous relations where the 
unknowns are the R-matrix entries a, b, b, c, d and f. We have a number of possibilities of 
selecting six out of nine equations to construct a consistent linear system for these variables. For 
a given system to have a non trivial solution the determinant of its coefficients depending on the 
weights of the Lax operators weights must vanish. We can for example choose Eqs. (I21H26P and the 
corresponding determinant can be written in the following form, 



/2 i n i ^ a Ji t^fi'Z ii ji N / ,/ 



(b — a j )c a — {b — a j )c a a a c j —babe A a b b — c bb =0. (30) 



/2, //— // //2, 



/ — / 
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From Eg. (1301) we see that in principle we have three possible branches to be analyzed depending 
on the factor we choose to vanish. However, by analyzing other possible systems of six equations 
we noted that the common factors shared by their determinants are only the first two terms of 
Eq.f l30p . As examples of alternative choices of systems we would like to mention those built up 
from either Eqs. fl2m25|29|) or Eqs. f l21tl24|26|29p . In any of theses cases the form of the third factor 
always changes and thus it plays the role of an extraneous term. From now on we shall disregard 
the last factor in Eq.(l30l) as feasible branch. 

We further notice that only the first factor of Eq. (!30|) has the suitable polynomial form (fllj) . In 
fact, the second term of Eq.f l30p clearly does not vanishes as w' — > w". This is not an impediment 
to define an analogous of a divisor but this will leads us to distinct components for each of the 
weights labels. At this point we can not discard the second factor as a possible branch since we 
have not yet studied the full properties of the linear system. A more detailed analysis shall reveal 
us that such apparent asymmetry of the second factor divisor disappears. 

4.2.1 Main Branch 

This branch is defined by imposing that the first polynomial factor of Eq. (130]) is zero. It follows 
that the corresponding divisor is, 

" 2 ^ ! =A, (31) 



cd 

We have now the necessary condition to solve Eqs.( li2"TH2T)|) by linear elimination of the variables 
a, b, b, c, d and f . At this stage it is sufficient to present their expressions in a nested form since 
we have not yet solved the full system of nine equations. By using Eq.( l3T|) for double primed labels 
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we find that these R-matrix entries can be written as, 

— = [b co a — a a c j ) I (Aic a c a 

ii a it i i a n i i n a 
bead — a a b c ) I (Aic a c a J, 



c 
b 



c 

b .Si /ii i ii . . . ' " \ 

— = [—0 c —cb)/[ac), 
c c 



d 

c 

f 

c 



7 / // , b / // b / // 

a c i—b a j 6 

c c 



/ 



iji n2 ii ii—n 

b b c — And b b 



b— / / ni b i2 n— a 
—b j c — An— d a b 
c c 



/ 



/-/ nl i2 iijn 

b b c — And b b 



(32) 
(33) 
(34) 
(35) 

(36) 



where c is an overall normalization. 

We see that Eqs (l32ll33|) are singular when the parameter Ai is zero and this means that here 
we have to assume Ai ^ 0. As we shall see the special case Ai = will be covered by the branch 
related to the second polynomial factor of the determinant ( l30l) . Let us now begin the analysis of 
the remaining three relations considering first Eq. (1271) . After substituting the expression for the 
ratio b/c (1331) in Eq. (l27|) we easily find that it becomes proportional to the polynomial, 



(be — Aib d)a d — (be — A\b d )a d = 

which again has the form (JTTJ) and the associated divisor is, 

be — Aibd 



ad 



A 



2- 



(37) 



(38) 



We next consider the solution of Eq.( l28p . We first observe that we have already been able to 
reduce the number of independent variables by two weights. In fact, considering Eqs.( l31|l38|) it is 
not difficult to resolve the weights / and d in terms of the remaining variables a, b,b and c. By using 
this information together with the expression for the ratio b/c (|34|) and after few simplifications 
we find that Eq.( l28|) can be expressed as follows, 



A 'V , A 'j/ 2 A ° V 2 

A 2 a o +Aiao — —ab 
Ax 

A " 2 U" i A "U' 2 ^° 'V 2 

A 2 a b +Aia b — —a b 

Ai 



. // , // 2 , , // 2—// . _ // // 2 

A 2 a b + Aib b - A x b c 



' i.i2 . ,2-i . , ,2 

A 2 ab + Aib b - A ± b c 



(39) 
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Notice that the polynomial (I39|) is already written in the form (TTTT) . Its dependence on the 
parameter Aq can be re-scaled by means of the following transformation, 



Ai = A1VA0 and A 2 = A 2 VA , (40) 

and now Eq. (l39l) becomes only dependent on Ai and A 2 . This re-scaling of parameters is ultimately 
due to the freedom of implementing a gauge transformation on the weights d and d. Taking into 
account Eq. fj40|) we find that the respective divisor associated to the polynomial fl39|) is given by, 



b 2 - AiA 2 a6 - Ajb 2 



A 2 ab + A1&6 - Aic 2 



A 3 . (41) 



We have now reached a point in which only Eq.( l29|) remains to be solved. As before we would 
like to use the last divisor ( 14T]) to eliminate one further weight and thus reducing the number 
of degrees of freedom. We observe however that the divisor (141 j) does not provide us the means 
to eliminate any weight in a linear way. This difficult can be circumvented once we analyze the 
polynomial associated to Eq.f l29|) and notice that it depends on the weight c only through even 
powers. This means that we can use the fact that divisor (|41|) has a quadratic dependence on 
the weight c to systematically eliminate this weight from Eq. (l29p . In Appendix A we explain 
how this operation can be implemented within the Mathematica algebraic computer system. As a 
result of this procedure we find an involved bihomogeneous polynomial with bidegree (4,4) which 
unfortunately can not be brought into the convenient form (11 II) . At this stage of the analysis the 
requirement that the polynomials must satisfy the anti-symmetric property ([TO]) becomes decisive 
to make further progress. By imposing this property we find that it can indeed be fulfilled provided 
that the so far free parameters are constrained by the following simple relation, 

A 3 A 2 - A? - 1 = 0. (42) 

After using the condition (142"j) we find that a large number of terms of the resulting polynomial 
coming from Eq. (l29"I) are magically canceled out. Thanks to this simplification we are able to 
write the resulting polynomial in the appropriate form (fTTj) . In what follows we shall present the 
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respective divisor since from it one can easily recover the associated polynomial. The expression 
for the divisor is given by, 



This completes the solution of the nine functional relations for this branch. Up to this point 
the effective integrable manifold of the weights should be given by the intersection of the last two 
divisors (j41|43p . In addition to that we have four free parameters at our disposal since so far we just 
have the constrain f )42|) . These divisors give rise to two projective surfaces G CP 3 and is generally 
expected that their intersection will lead us to an algebraic spatial curve. Since some of the free 
parameters are going to be fixed later on we shall postpone the analysis of the intersection until 
the very end. 

4.2.2 Special Branch 

This branch is defined by setting the second term of Eq.(l30l) to zero and it is related to the 
particular value Ai = excluded in the main branch. This can be seen by first noticing that the 
R-matrix ratio a/c ( 132]) is proportional to the second factor of the determinant (130]) . For this ratio 
to be no null we should have another zero on the denominator such that the value of the ratio 
becomes indeterminate. Inspecting Eq.( l32j) we note that the only option is to set Ai = and as 
result it follows that the weight / must be fixed by the expression, 



We now observe that the numerator of Eq. (133]) must vanish otherwise we would have a divergence 
on the ratio b/c. From this condition it follows a polynomial relation of the form (ITT]) , namely 



A 2 a 2 (A 2 b 2 - A 2 b 2 ) + A 1 A 2 b 3 (2A 2 a + Ai&) + b 2 (A 3 3 b 2 - A 2 b 2 ) 

6& 2 (Ai& + A 2 a) 



A 4 . 



(43) 




(44) 



I l ll ll ll ll l J 

bead — o c a a = 0, 



(45) 



whose corresponding divisor is, 




(46) 
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The next step is to evaluate the indeterminacy of the above mentioned R-matrix ratios. This 
is fortunately done with the help of Eqs. (j21|27p and the final expressions for such ratios are, 

a/c = [b a b — abb + a c \/{c c a ), (47) 

— , ,— I II I— II , / I If \ / V 

b/c= (ba -ab )/(cc ), (48) 



while the other R-matrix entries can again be computed from Eqs. (I34H 36]) . 

At this stage the only relations that remain to be solved are Eqs. (1281129]) . The technical details 
entering their solution are fairly parallel to those already explained in the previous subsection for the 
main branch. In what follows we shall therefore present only the final results for the corresponding 
divisors. We find that the divisor associated to Eq.( l28]) has the following form, 

A 2 a 2 b , , , 

rai = (49 > 

where A2 has been re-scaled as in Eq. (T4T)]) and together with the parameter A3 they satisfy the 
constraint, 

A3A2 -1 = 0. (50) 
Finally, the divisor associated to the anti-symmetrical polynomial derived from Eq. (129]) is given 

by, 

A 2 a 2 (A 2 b 2 - Ajb 2 ) + b 2 (A 3 3 b 2 - A 2 b 2 ) = A 

A 2 bb 2 4 ' 1 } 

We note that the divisors (I44f46f5ip can be directly obtained from those derived for the main 

branch by substituting Ai = in Eqs. (l31|l38P3]) . respectively. The same observation for the re- 
maining divisor (149 p is more subtle since the zero order in Ai of the corresponding main branch 
divisor ( 14"T]) is trivial. In spite of that we can recover the divisor (]4"9~]) by rewriting Eq. (l4ip in powers 
of the parameter Ai. By setting the coefficient proportional to Ai to zero we then easily obtain the 
divisor ( T41"]) . In practice, since this limit is somewhat delicate we shall consider these two branches 
separately. 
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4.3 Four Terms Equations 

The Boltzmann weights h, h and g begin to emerge only in the functional equations having 
four distinct terms. Even after using the divisor ( 12 Op we find that the total number of relations 
with four terms is still considerable. Altogether we have twenty-two functional equations which 
is a high number to approach the problem by the standard elimination method. As we shall see 
however we need to solve only ten functional relations to decide on the integrability of the vertex 
model. The reason for this simplification is that few of the equations with four terms have as 
unknowns the R-matrix entries a, b and b which have already been fixed. For sake of consistent of 
the elimination procedure we are then forced to make linear combinations between a small subset 
of four terms functional equations and two particular three terms relations. Remarkably enough 
this simple reasoning is able to determine the structure of the remaining weights h, h and g. We 
shall first detail how this procedure works for the weight g. 

4.3.1 The weight g 

We have three functional equations depending uniquely on the weight g and on some variables 
that have been previously determined. Their expressions are given by, 

eg o + be c — bet a — co g = 0, (52) 
get b + do c — co a —hag = 0, (53) 
go c + dab — cc o — hg c =0. (54) 

We see that Eqs.( l52ll54l) have five unknowns and therefore we need more two equations to build 
up a consistent homogeneous linear system. These extra relations should be searched among those 
solved in subsection (14.2. ip since the above equations depend on the weights b, b and c. Direct 
inspection of Eqs.( l21t l 2^|) reveals us that this choice is remarkably unique once we want to keep the 
minimal number of five unknowns. These suitable relations turn out to be Eqs.( l27|l28|) . Now, by 
setting the determinant of equations (I27ll28|52tl54|) equal to zero we find that it can be factorized 
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as, 



y_/ / 1 / // 1 / / // _ I I II 1 f f I / v k ' J / T ' l' ^ " T" ^ T ' ' ' t'" 

A (o c a c a g — bcdcdg) + A (o dad — b d ad 

2, //-//" 



-/ /2 //2 // -// //2 /2 
+OCC / — OC C J 



u'u' " 2 
b b c 



And b b 



0. 



(55) 



The determinant ( 1551) must vanish through the first factor since the second one is exactly the 
extraneous term that has been discarded before. We see that the first factor of Eq. (1551) is anti- 
symmetrical on the exchange of weights labels and taking into account the divisors (I31|38|4ip we 
are indeed able to rewrite it in the appropriate form ( TIT]) . The resulting polynomial depends only 
on the weights a, b ,b, g and the corresponding divisor is given by, 

A 2 A 3 6(A 2 a + Aib)g - b[A 3 b 2 + A 2 (A 2 a + A^) 21 



6 2 (A 2 a + A 1 6) 



A 5 



(56) 



Let us consider the solution of the remaining equations ( !53ll54l) . We can easily solve one of these 
relations fixing the value of the R-matrix element g. For instance from Eq.( l53|) we find, 



b / // — / // d / n 

—d g + b d be 

c c 



/(d'b" 



(57) 



where the ratios b/c and d/c have already been determined by Eqs.( |34ll35i) . 

The last equation ( 15^1) can now be rewritten only in terms of the weights a, b and b once we 
use the help of the divisors ( |3T|38II41|[56|) . After some simplifications we find that Eq. (|54j) becomes 
proportional to the following expression, 



1 2 



A 5 - Ai(Aa + A 2 ) I l6"'(A 2 A 2 a " + A 2 6'" + A 2 A 5 a b' + (A^ - A\)b' 2 )+ 
A 2 (A 3 - A 2 )a"b"(A 2 a'b' + A^' 2 ) - A 2 6 7,2 (A 2 A 3 a' 2 + Ai(A 3 + A 2 )a'&' + A 2 6' 2 ) 



(58) 

The polynomial (1581) is very far from satisfying the anti-symmetrical property (TlQi) but never- 
theless it is proportional to a combination of free parameters. We then are able to solve Eq.([5i 
by choosing that its first factor vanishes and as result the parameter A 5 becomes fixed by, 



A 5 = A!(A 3 + A 2 ) 



(59) 



We finally remark that though the above discussion has been detailed for main branch there is 
no difficulty to repeat the same analysis in the case of the special branch. It turns out that for the 
special branch we simply have to set A x = in the final results 
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4.3.2 Weights h and h 

Out of thirteen functional equations involving the weights h and h we shall need only seven of 
them to determine both these weights and the R-matrix elements h and h. The expressions of such 
basic relations are, 



, Ijll f- / II 

chb + be c 


— hdd — cb h = 


0. 


(60) 


7 ij" f ' " 

ch o + be c 


— hdd — cb h = 


o, 


(61) 


hb c — hh c 


— ccb + da o = 


0. 


(62) 


hb c — bh c 


'," ^ijii 
— ccb + da o = 


0. 


(63) 


ha b + do c 


— cb d —hdh = 


o, 


(64) 


ha b + do c 


— cb d —hdh = 


o, 


(65) 


ah a — dc d 


— fh f — ha h - 


= 0. 


(66) 



As before we observe that Eqs. (160116 ip have as unknowns the R-matrix elements a, b and b and 
these equations can again be solved by making linear combinations with Eqs. (127|28|) . By requiring 
that the determinant of Eqs. fl27|28|l60|) vanishes we obtain, 



cd c d 



b'(h" -a" -f/A )-b"(ti -a -f/Ao) 



0. 



(67) 



while the vanishing of the one made by system of equations ( J27|28f6ip is, 



c dc d 



b'(h" - a" - A f") - b"(fi - a - A Q f) 



a . — // . — / / 



0. 



(68) 



The second factors of the above relations give rise to anti-symmetrical polynomials of the form 
( II ip and their corresponding divisors determine the expressions for the weights h and h, 



h-a-f/A h-a-Aof x 
= = A fi and = = A f 



(69) 



We next note that from Eqs. 



we can retrieve the R-matrix entries h and h. We further 



observe that the last two terms of these equations are also present in the previously solved three 
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terms relations (I22|I26[) . This makes it possible to eliminate the terms cc b and dd b of Eqs. (162(63]) 
and after using the divisors (|69|) for the single primed labels we find that the ratios h/c h/c are 
given by, 



a .f 6 b 

- + (- /A + A 6 -- 
c c be 



and 



a A /f\ 

- + A o 
c c 



A 6b 

A 6 77- 

c 



(70) 



h 

C C C DC c 

We now have the basic ingredients to tackle the solution of Eqs. ( 16411651) . Considering the above 
results altogether as well as the previous expressions for the ratios d/c and b/c we are able to write 
these equations solely in terms of the weights a, b and b. We find that the vanishing of Eqs. (I64f65p 
are equivalent to the the following identities, 



A 6 A &' A 2 (b' 2 a"b" -b""a'b')+A 1 (b"'b""-b'"'b") + A 6 A 3 6'"o" (a b" - a" b' ) = 0, (71) 

A 6 V [A 2 (6' 2 aT - b" 2 a'V) + A^&'V 2 - b" 2 b' 2 )] + A 6 A 3 A 6'Y (a'b" - a'T) = 0. (72) 
We see that Eq. (l71|72l) contain common anti-symmetrical polynomials which in principle could 



b a b 



>2 b" 2 -b" 2 b' 2 ) 



2—// . i— ii 



be set to zero by means of the corresponding divisors. This possibility imposes at least an extra 
constrain on the weights a and b in addition to the other two divisors ( 1411431) already derived in 
subsection ( 14.2. ip . The intersection of such three divisors will generically lead us to zero dimensional 
manifold consisted of finite number of points for the ratios a/c, b/c and b/c and thus to Lax 
operators without free spectral parameters. From now on we shall discard this kind of possible 
"braid" solutions of the Yang-Baxter algebra. However, it is fortunate that Eqs. (171l72j) can vanish 
without the definition of any additional divisor by choosing the last two free parameters to be zero, 
that is, 

A 6 = A 6 = (73) 

Let us finally consider the solution of Eq.( l66p . Considering all the results we have obtained so 
far we find that Eq.( l66p leads us once again to deal with a very complicated polynomial of bidegree 
(4,4) on the weights a, b and b. To make progress we then impose the anti-symmetrical property 
(flUI) in the expectation of further simplifications. Fortunately, we find that such property can be 
satisfied provided we introduce a new constrain among the parameters, namely 



A 4 AiA + A3 



A 2 (1 + A2)-A 3 A ) - A 2 A (1 - A?) = 0. 



(74) 
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By substituting the constrain (17^1) back into Eq. (!66|) we find that its non trivial terms become 
proportional to the following polynomial, 



i— n a— >i. 



ii t til ,~//2v — // / r ul 



(Ao-Ao+l)(a b -a b ) A x a (b + b' ) + b (A 2 a + A 3 6 ) A x a (6 + 6" ) + 6 (A 2 a + A 3 6 

(75) 

Taking into account the above discussion we conclude that Eq. (I66[) should vanishes by imposing 
that the first factor of (175]) is zero. As a consequence the parameter Aq is constrained to take values 
on the non trivial third order roots of unity, namely 



7T, 



A = exp(±i-). 



(76) 



Now by considering the above constrain for Ao back to Eq.( IT4l) we see that this parameter 
factorizes leading us to an effective condition independent of A , 



A 4 A X + A^(A 2 - A 3 ) - A 2 (l - A?) = 0. 



(77) 



Although the above discussion has been concentrated in the case of the main branch the same 
reasoning applies straightforwardly for the special branch. Once again we just have to impose 
Ai = in the final results and we see that in practice this becomes relevant only for Eq. (1771) . By 
substituting Ai = in this relation and with the help of Eq. (l50p we conclude that for the special 
branch the parameter A 2 is required to satisfy the relation A 2 — A| + 1 = 0. This means that the 
for special branch all the parameters have already been fixed with the exception of A4. By way of 
contrast for the main branch we have two free parameters and to avoid division by a potential zero 
factor we choose them to be A 2 and A 3 . For sake of clearness we have summarized the values of 
fixed parameters for both branches in Table [2 



Parameters 


Ao 


Ai 


A 2 


A3 


A 4 


A 5 


A 6 ,A 6 


Main Branch 


exp(±i|) 


i^A.As - 1 


Free 


Free 


A|(A 3 -A2)+A 2 (1-A2) 


Ai(A 3 + A 2 ) 





At 


Special Branch 


exp(±i|) 





+1 
+A 2 


A 2 1 


Free 









Table 2: The values of the parameters associated to the divisors for both branches. 
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We now reach a point in which there are twenty-four functional relations involving four and five 
terms that still have to be verified. To avoid overcrowding this section with a number of additional 
formulae we have presented their explicit expressions in Appendix A. Taking into account the 
results obtained so far the bihomogeneous polynomials associated to Eqs. (IA.i]|A.24j) can clearly be 
expressed in terms of the weights a, b, b and c. In the case of the main branch we can use the 
help of the divisors f j4T|43|) to simplify these polynomials in a systematic way until we reach the 
point where the weight c is completely eliminated and the powers of the weights b are at most 
three. Remarkably enough, after such algebraic manipulations these twenty-four polynomials are 
either zero or become proportional to the factor Aq — Ao + 1 and consequently vanish thanks to the 
constrain (1761) . Similar reasoning can be implemented for the special branch by considering now the 
respective divisors (I49|51l) . The technical details concerning these simplifications are described in 
Appendix A including the explicit computer algebra system code used. In this way we are able to 
verify algebraically that the whole Yang-Baxter algebra (J5j) is indeed satisfied for both branches. 
We emphasize that this verification can be done in rather modest computer as far as memory is 
concerned and without the need to take any a priori numeric values either for the available free 
parameters or for the Boltzmann weights. 



5 The Manifolds Geometry 

The purpose of this section is to present the structure of the main algebraic manifolds governing 
the integrability of the previously uncovered vertex models as well as to discuss their geometric 
properties. We shall argue that the weights of the Lax operators can be written solely in terms of 
three variables a, b and c which are constrained by projective plane algebraic curves of genus five 
for generic values of the free parameters. Since we will be dealing mostly with singular manifolds 
we start by recalling the definition of singular locus of an irreducible hypersurface S(u , . . . , u m ) G 
CP m . The set of singular points of this hypersurface form a closed subvariety Sing(S') determined 
by the zeroes of all the partial derivatives of S(cuq, . . . , cu m ), 

Sing(S) = {[u :...:u> m ]e CF m \— = 0, for ,j = 0, . . . ,m.} (78) 

duij 
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The analysis of the singularities is going to be helpful to elucidate the intersection problem of 
surfaces in the case of the main branch. This study is also essential for the understanding of the 
geometric properties of the algebraic curves we shall confront here. 

5.1 The Main Branch 

From the majority of the divisors we are able to extract the weights in terms of the variables 
a, b and c by means of linear elimination. The only exception for this branch is the weight b 
since we have to deal with the intersection of the divisors f)41|43p which are both non-linear in this 
variable. Before considering this problem we find convenient to perform the following re-scaling on 
the parameters A 2 and A 3 , 

A 2 = Ai£i and A 3 = A 1 e 2 - (79) 

where e x and e 2 are now the two free parameters associated to the main branch. 

Considering the above re-scaling and the parameters values of Table ([2]) we can rewrite the 
divisors f)41|43p as two projective surfaces in CP 3 . The polynomial expression associated to divisor 
(j4"Tl) becomes, 

Si(a,6,6,c) = a(b 2 + b 2 ) + (e 1 a 2 + e 2 b 2 )b-e 2 bc 2 , (80) 

while the divisor ( l4"3l can be rewritten as, 

S 2 (a, b, b) = a 2 {elb 2 - e\b 2 ) + (2e x a + b)b 3 - {e x e 2 -e\- 2)b 2 b 2 - {e 1 e 2 - l)b 4 

- ie\e 2 - e x e\ - 2e x + e 3 2 )ab 2 b. (81) 

The surface Si is free of singularities for generic values of the parameters e\ and e 2 because there 
is no solution to Eq. (178p other than the origin. This is a rational manifold since any nonsingular 
cubic surface is known to be birational to CP 2 [H]. The surface S 2 does not depend on the weight 
c and thus has to be seen as cone over a plane algebraic curve. By solving Eq. (1781) for S 2 we find 
that it contains two non coplanar singular lines given by, 

li = [a : : : 1] and 1 2 = [a : : -e x a : 1], (82) 
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and for this reason it is expected that S 2 should be a cone over an elliptic curve [T7]. In appendix 
B we present the technical details showing that £2 indeed defines a curve of genus one when viewed 
in the ring C[o, b, b]. 

Let us turn to the problem of computing the intersection of the pair of algebraic surfaces (180118 II) . 
From the algebraic geometry theory we know that the result of the intersection of two irreducible 
surfaces Si and S2 C CP 3 is always a finite collection of curves [18] . An extension of the Bezout's 
theorem asserts that the product of the degrees of Si and S 2 should satisfy the relation, 

deg{S l )deg(S 2 ) = Y,USi,S 2 )C a . (83) 

a 

where the sum varies over all the irreducible curves C a of Si fl S 2 . The index l a (Si, S 2 ) is the 
multiplicity of the intersection of Si and S 2 along the component C Q . 

The first step of the problem is therefore to find out whether the intersection is irreducible or 
it degenerates in the union of a given number of curves. This can be answered by exploring a basic 
fact in commutative algebra that assures us that the ideal generated by the polynomials Si and 
S 2 can be decomposed in terms of more elementary ideal components, see [19] for mathematical 
details. These components are called primary ideals whose defining polynomials give rise to the 
varieties corresponding precisely to the irreducible curves coming from the intersection Si H S 2 - It 
is fortunate that the existing algorithms for extracting information on this primary decomposition 
have already been implemented in computer algebra systems such as for example Singular [20J. 
With the help of this software we find that the ideal associated to the polynomials (I80II81I) is in 
fact reducible having three main primary components. Here we are in a favorable situation in 
which two of them sit exactly at the singular locus of S2 defined by the lines (1821) . By analyzing 
the factorization of the polynomials (I80|l8ip around b = we can easily obtain that the respective 
intersection indices are Ii^SV, S 2 ) = Ii 2 (Si, S 2 ) = 2. Considering all these information together with 
formula (1531) we conclude that the third component has to be a spatial curve of degree eight. From 
the practical point of view such third component is the only non trivial result of the intersection 
since we are not interested in solutions having null weights. 

The next step is to map the surface intersection to a degree eight plane curve by hopefully 
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determining exactly one of the weights in terms of the curve variables. To this end the result of 
the primary decomposition for the third component is of any help since the respective basis of 
the primary ideal is constructed out of a large number of polynomials. This task can however be 
resolved by means of the following method. We rewrite the the surfaces expressions as univariate 
polynomials in the variable we want to eliminate which here we choose to be the weight b, namely 

Si(b) = v 2 6 2 + Vl 6 + v , (84) 
S 2 (b) = u 4 6 4 + u 3 6 3 + u 2 b 2 + U16 + u , (85) 

where the coefficients Vj and Uj depend on the remaining weights a, b and c. They are determined 
by just matching Eqs.( 180f8T|l to Eqs. (l84ll85l) respectively. 

We now proceed by making certain linear combinations among S\(b) and 82(b) in order to lower 
the highest degree in the variable b. For instance this is achieved by defining the following new 
polynomials, 

Sx(b) = [u St(b)-v S 2 {b)]/b, 

S 2 (b) = u 4 b 2 S 1 (b)-v 2 S 2 (b). (86) 

Because both Si(b) and S 2 (b) are algebraic combinations of the starting surfaces they certainly 
contain the intersection curve that we are searching for. The nice feature of the combination ( |86j) 
is that the maximal degree of the polynomials has now decreased to three. By repeating this 
procedure using now in the right hand side of Eq.f l86|) the new polynomials Si(t>) and S 2 (b) we are 
able to lower the degree in b once again by one. After three such steps we reach a couple of relations 
that are both linear in b and this variable can finally be eliminated. By extracting the weight b 
from one of the relations and substituting the result back into the other we obtain a polynomial 
in the variables a, b and c that factorizes in terms of three curves with degrees six, eight and nine. 
The components with degree six and nine play the role of extraneous factors and they should be 
discarded since our previous rigorous argument tell us that the sought plane curve should have 
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degree eight. The explicit expression for this algebraic plane octic curve is, 

Ca(a, b, c) = (eie 2 - l)[a 4 + a 2 b 2 + b 4 } 2 + {2 + e\- e x e 2 + e%) [e 2 {a 4 + a 2 b 2 + 6 4 ) + abc 2 } abc 2 
- [(e x e 2 - 2)a 4 + (2 - e x e 2 + £ 2 )& 4 - £ 2 a 2 fc 2 + 2e 2 abc 2 + c 4 ] c 4 , (87) 

and after a systematic use of the constrain (IBT1) we are able to simplify the expression for the weight 
b to obtain, 

- = [(e 2 - gi)Q 2 + g 2 ^ 2 ] be 2 - a[a 4 + a 2 fr 2 + b 4 - c 4 ] 

e 2 (a 4 + a 2 b 2 + 6 4 ) + 2abc 2 ' 1 j 

As possible check of the above results we can substitute the weight (I88p in the original polynomial 
equations defining the surfaces Si and S 2 . With the help of a symbolic algebra system one easily 
find that the relations (180M81I) become indeed proportional to the octic plane curve ( |87l) . 

We now turn our attention to discuss the geometric properties of the degree eight plane curve 
(IB71) governing the integrability of the main branch manifold. The basic invariant characterizing a 
projective algebraic plane curve is the topological genus of the respective compact Riemann surface 
normalization. In order to compute the genus we need first to identify the singular points and 
afterwards to investigate their morphology which include the understanding of possible infinitesimal 
neighboring singularities. By solving the polynomial equations f JTHj) for the plane curve flHTl) we 
find a total number of twelve singular points. The first eight of them are located in the affine plane 
c = 1 and they can be expressed as, 

P A =[a a :6 a :l], (89) 
where the values a s and b s are the non-null solutions of the following relations, 

e 2 2 al + (ei£ 2 - 2)a s b 2 + e 2 b s = 0, (90) 
e\bl + (2 - e x e 2 + e 2 2 )a 2 b s + e 2 a s = 0. (91) 

The above coupled non-linear equations can be resolved in terms of the roots of a single uni- 
variate polynomial with degree eight. The final answer is somewhat cumbersome but for sake of 
completeness it has been presented in Appendix C. At this point we recall that a relevant index 
characterizing the morphology a singular point is the multiplicity of the singularity. A singular 
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point [ao : &o : cq] is said to have multiplicity m on the plane curve Ci(a, 6, c) if every partial deriva- 
tive with respect to the coordinates a, b and c up to the order m — 1 vanishes at [do : &o : Co] while 
at least one of order m is non null at this point. In addition, the singular point is called ordinary 
when the factorization of the leading term of the expansion of Ci (a, b, c) around [a : b : cq] gives 
rise to m different terms and thus we dont have multiple tangents through this point. It turns 
out that the Taylor expansion of Ci(a, b, 1) near to the affine singularities defined by Eqs. ( I891I9T1) 
always produces two distinct tangents for generic values of the parameters e\ and 6%. This means 
that all the affine singularities are in fact ordinary double singular points. 

The remaining four singular points are located at the infinity line c = being just the zeros of 
the monomial C\(a,b, 0). As a consequence they are independent of E\ and e<i and their explicit 
coordinates are, 



The singularities at infinity behave as double point with only one tangent but having two 
branches. They are not ordinary singular points and are usually known in the literature as tac- 
nodes. Their presence means the existence of extra neighboring singular points and the plane 
curve Ci(a, b, c) desingularizes only after the implementation of a sequence of two global birational 
transformations often named blowing-ups. The main idea of this method is to replace a point in 
the plane by a projective line which opens more room for the curve to become non-singular in 
a higher dimensional space. For the technical details concerning this approach in the context of 
desingularization of an algebraic plane curve we refer to an overview by Abhyankar [21]. As a result 
the resolution of the singularities of the octic plane curve (j57j) can be represented by the diagram, 



where 7i"i and 7^2 represent two consecutive blowing-ups, the intermediate curve Ci carries infinitely 
near singularities and the smooth normalization is denoted by Ci. 

Under the transformation 71*1 all the affine singular points are mapped onto simple non-singular 
points while the singularities at the infinity line become four ordinary double points. The latter 
are called infinitesimal neighboring singular points now sited in C\ which are finally resolved with 




and 




(92) 




(93) 
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the help of the second map 1^2- The genus of the curve normalization can be computed using the 
following standard formula valid for plane curves, 

9(C0 = (D-mo- 2) (94) 

p 

where D is the degree of the curve and mp denotes the multiplicity of the singular point P. The 
sum is taken over all the singular points on the curves that are infinitesimal neighbors of C\(a, b, c). 

At this point we know that all the singular points have multiplicity mp = 2 including those in the 
infinitesimal neighborhood of the singular points at the infinity line. Considering this information 
in Eq.f l94|) we obtain that the genus is, 

g{Cx) = ^-12^L-4x_L^=5, (95) 

7T1 7T2 

whose value can be confirmed within symbolic algebra packages capable of computing the genus of 
plane curves such as Singular [20J. 

In general, canonical genus five curves are known to be realized as the complete intersection of 
three quadrics in CP 4 but they also can degenerate in either hyperelliptic or trigonal curves, see for 
example [22]. One effective way to shed some light on the actual class of the algebraic curve ( 1871) 
is to investigate possible mappings to other curves with lower genus. Exploring the fact that the 
polynomial ( 187)) depends only on even powers of the weight c we are able to establish the following 
regular map, 

Ci(o, b, c) C CP 2 -A Qi(x, y, z) C CP 2 

(96) 

[a : b : c] 1 — > [a 2 : ab : c 2 ] , 

where the image of the map <ft is the algebraic curve Qi(x, y, z) defined by, 

Qi(x, y, z) = (6x62 - l)[x 4 + x 2 y 2 + y 4 } 2 + (2 + e\ - e x e 2 + e 2 2 )[e 2 (x 4 + x 2 y 2 + y 4 ) + x 2 yz]x 2 yz 
- [(s 1 s 2 -2)x 4 + (2-e 1 e2 + e 2 2 )y 4 -e 2 2 x 2 y 2 + 2e2X 2 yz + x 2 z 2 ]x 2 z 2 . (97) 

The degree of the above map is two because this is the cardinality of the fiber _1 (P) for a 
generic point P G Qi(x, y,z). To make further progress on this admissible double cover we need 

29 



to compute the genus of the target curve. In this sense we find that Qi(x, y,z) has eight ordinary- 
singular points and one singularity resembling the tacnode behaviour but with higher multiplicity 
mp = 4. For sake of completeness the technical details entering this analysis have been summarized 
in Appendix C. The curve Qi(x, y, z) desingularizes once again after a sequence of two blowing-ups 
and the desingularization diagram is similar to that shown in (193]) . Denoting the corresponding 
blowing-ups by 7Ti and ^2 we find that genus of this curve is given by, 



where Qi denotes the normalization of Qi(x, y, z) which turns out to be an elliptic curve. 

Now, putting all these information together we can establish the following commutative diagram, 



where ifj denotes the morphism induced by the map 0. 

From the diagram (199 p we see that the morphism ip has also degree two and therefore the genus 
five curve Ci admits a double covering in terms of the smooth elliptic curve Qi. The existence 
of such degree two map is known to rule out the possibility that the octic plane curve fl87|) be 
sited in the space of either hyperelliptic or trigonal curves [23] ■ We then conclude that we are in 
fact dealing with a bielliptic genus five curve whose canonical model is the complete intersection of 
three independent quadrics in CP 4 . 

We close this subsection emphasizing that the above discussion on the geometric properties of 
the octic plane curve ( 18 7p is valid for generic points on the two-dimensional space generated by 
the free parameters E\ and 62- Such family of curves may degenerated to plane curves with lower 
genus once we restrict S\ and 62 to lie on certain specific submanifolds of the parameter space. This 
fact can naturally occur when one of the surfaces (I80|I81I) becomes reducible and as result their 
intersection will give rise to reducible plane curves whose components should have degree lower 




(98) 




(99) 
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than eight. In Appendix B we have in fact remarked that the the surface (IHTj) can degenerate 
into the product of two cones over plane conies. This happens when the parameters e\ and £ 2 are 
constrained to sit in the following one-dimensional submanifolds, 

e 2 exp(±i^)-e 1 exp( T i^) T 2i = (100) 

O D 

£ 2 exp(±iJ)- £l exp( T iJ)±2i = (101) 
o o 

Ae\ + e\- 2e\e 2 + te\ + Ze\e\ - 2e x e^ + e 4 2 = 0. (102) 

It is not difficult to repeat the previous reasoning on the intersection of the surfaces (I80|l8ip 
now for the particular submanifolds (I100H1021) . Due the factorization of the degree four surface (I8ip 
one expects that the result of the intersection will lead us to the product of two quartic algebraic 
plane curves. This analysis is somehow direct for the first two manifolds ( llOOOlOlj) since one of the 
parameter can be linearly eliminated but it is more involved for the third non-linear submanifold 
(11021) . The technical details entering this analysis are presented in Appendix D and here we only 
state our main conclusions. From the explicit expressions of these quartic plane curves we conclude 
that they are singular. The singularities associated to the linear submanifolds ( II 0011 10 ip are single 
tacnodes at the infinity line while those related to the third submanifold ( I102p are constituted of 
two ordinary double points in the affine plane. This means that when the parameters E\ and e 2 are 
constrained in any of the special submanifolds ( 1100111021) the main branch integrable vertex model 
is therefore governed by plane quartic curves of genus one. 

5.2 The Special Branch 

The situation for the special branch is much simpler since the weight b can be linearly eliminated 
from the divisor (H9|) . Considering the results of Table ([2]) one would think that this branch has 
to be splited in four different integrable manifolds associated to the four possible values of the 
parameter A 2 = ±exp(±i|). However, these apparent distinct vertex models are related to each 
other and as result we have only one independent integrable manifold. The multiplicative signs 
are easily gauged away with the help of standard gauge transformations. In addition, the vertex 
models with A 2 = exp(i|) and A 2 = exp(— i|) can then be connected by applying the Weyl basis 
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transformation 1 f)- 3 on the Lax operator (1171) . Taking this observation into account we find that 
the expression of the weight b for such single vertex model is, 

An&c 2 



a 2 + A 6 2: 

where we recall that Aq = exp(±i : |). 



(103) 



In order to obtain the constrain among the variables a, b and c we substitute the above expression 
for the weight b in the divisor (l5"Tj) . After some simplifications using the identities satisfied by the 
parameter A we find the respective algebraic plane curve is given by, 

C 2 (a, b, c) = [a 2 + A b 2 ](a 4 + a 2 b 2 + A A abc 2 ) + {b 2 - a 2 )c 4 , (104) 

where A4 is the free parameter. 

The sextic plane curve (I104p has three singular points being one of them on the affine plane 
while the others are sitting on the line at infinity. Their explicit coordinates are given by, 

P A = [0 : : 1], Poo = [i- : 1 : 0] and [~ : 1 : 0]. (105) 

Ao Ao 

The singularity at the origin of the affine plane is an ordinary double point while the remaining 
ones behave as tacnodes and the sextic curve is desingularized again by a sequence of two blowing- 
ups. We can compute the genus of the normalization C2 of the sextic curve 02(0, b, c) along the 
same lines presented in previous subsection. By applying the formula (j94p we obtain that the genus 
of C 2 is, 

g (C 2 ) = - 3_x_l-2_x_l = 5. (106) 

TTl 7T2 

In what follows we shall argue that the sextic plane curve turns out to be a bielliptic genus five 
curve as well. The first step is to note that the same two sheeted cover we have discussed before 
can also be established for the sextic curve ( I104p . namely 

C 2 (a, b, c) c CP 2 A Q 2 (x, y, z) c CP 2 

(107) 

[a : b : c] 1 — > [a 2 : ab : c 2 ] , 
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where the expression of the algebraic target curve Q 2 (x,y, z) is, 



Q 2 (>, y, z) = [x 2 + A y 2 ](x A + x 2 y 2 + y A + A 4 x 2 yz) + (y 2 - x 2 )x 2 z 2 



108) 



We next observe that the image curve (11081) is quadratic in the variable z and therefore the 
linear term on this variable can be eliminated by quadrature. In an analogy to what has been 
explained in Appendix B this define a birational transformation now in the projective space. More 
precisely, we are able to put forward the following second mapping, 



Q 2 (x,y,z) cCF 2 
[x:y:z] 



Q 2 (a?i, 2/1,21) C 



4>i(x, y, z) : fa(x, y, z) : 4> 3 (x, y, z) 



where the polynomial components of the map H109[) are given by, 



bx(x,y,z) = x(x +A y ), 
l 2 (x,y,z) = y(x 2 + A y 2 ), 
h(x,y,z) ■ 



-i [(A 4 y(x 2 + A y 2 ) + 2(y 2 -x 2 )z] 



while the expression of the image curve Q 2 (xi, yi, z\) is, 



_ 4 

Q 2 {x 1 ,y 1 , zi) = x\z\ - —y\ - (4A - A 2 ^x\y\ + ±x\ 

^0 



(109) 



110) 



111) 



Direct inspection of the quartic plane curve fillip reveals us that it can readily be brought into 
the Jacobi's form of an elliptic curve. Considering the combination of the above two maps we are 
able to build up the following diagram, 



C 2 (a,6,c) 



Q,2(x, y,z) 



\ 



112) 



Q 2 (x 1 ,y 1 ,z 1 ) 



The degree of the composition of transformations among varieties having same dimension is the 
product of the degrees of the individual mappings and from this fact it follows that deg(0 o(p) — 2. 
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This means that the diagram f l 1 1 2 j) represents a direct double cover mapping from the singular 
genus five curve 02(0,6, c) to a quartic singular curve with genus one. As before this shows that 
the sextic plane curve (I104p sits in the moduli space of bielliptic curves. 

Finally, we remark that another direct consequence of the diagram (11121) is that at the special 
values A4 = 12A , — 4A the sextic curve (I104p should degenerate into a plane curve of genus one. 
This is because at these parameter values the discriminant of the elliptic curve (^2(^1, yi, z\) is zero 
and as result we have a double cover mapping 0o0 from the original plane sextic curve to a rational 
conic curve. In fact, for A| = 12A and — 4A we verified that sextic curve 02(0,6, c) factorizes in 
terms of the product of two non-singular cubic plane curves. 

6 The Integrable Lattice Models 

We now gathered the basic ingredients to present our main results in terms of the language often 
used in the modern algebraic theory of exactly solvable lattice systems. At this point we know that 
the Boltzmann weights of the uncovered vertex models sit on algebraic plane curves and in principle 
they can be uniformized by means of a single spectral parameter. For future convenience we rewrite 
the Yang-Baxter algebra in the following more general form, 

R a>/3 (A, /i)L Qifc (A)L /3>fc (/i) = L /9j& (/^)L Q!)fc (A)R Q!)/ g(A, fi), k = 1, • • • , N, (113) 

where a, (3 indicate the R-matrix action on the auxiliary spaces while the index k represents the 
quantum spaces. The spectral parameters are denoted by A and \x. 

We start by presenting the explicit expressions for the Lax operators L a fc (A) and the R-matrix 

Rq,/3(A, /i). 
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6.1 Lax operator and R-matrix 

Without loss of generality we shall normalize the Lax operator by the weight c. It turns out 
that for both branches the Lax operator (TTTj) can be written as follows, 



L a ,£;(A) 



a(A)[e£> ® ei? + <&> ® e&'] + 6(A) [eff 3 gW + f 

,(«) I „(«) O I „(«) <» I J Q ) O „(*) 



» 



,(«) 



,(*) 



(a) 



3 Wi 

'22 J 



6(A)[e£ ) ®ei K 1 



'33 



'22 



® e 



(*)i 

33 J 



(a) 



9 (*0 
'32 



+ e 



(«) 

23 



e£ } ] 



+ exp(±i-)d(A)[e^ ) ® eg + e£ } ® eg] + /(A)[eg> ® eg + eg ® eg] + 0(A) [eg ® eg] 



+ [a(A)+exp(Ti-)/(A)][e 13 



eg] + [a(A)+exp(±4)/(A)][e£ ) ®e {fc " 



13 
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The structure of the weights 6(A), d(X), /(A) and g(A) are however branch dependent as well 
as the underlying algebraic plane curves. For the main branch the weight 6(A) has been already 
determined by Eq.( l88l) while the remaining weights are obtained by solving the divisors (I31f38f56p . 
Considering the constrains among the parameters listed in Table fl2]) we find that the final results 
for these weights are: 
• The main branch 

[(e 2 - £i)a 2 (A) + £2 6 2 (A)] 6(A) - [a 4 (A) + a 2 (A)6 2 (A) + 6 4 (A) - l]a(A) 



6(A) 



d(A) 



g(x) 



£ 2 [a 4 (A) + a 2 (A)6 2 (A) + 6 4 (A)] + 2a(A)6(A) 



(115) 



eie 2 - 1 



6(A) 



exp(±if) [ ei o(A) + 6(A)]' 



/(A) = 6(A)[ £l a(A)6(A) + 6(A)6(A)-l] ^ 



a{X){e l a{X)+b{\)} 



eis 2 



6(A) + 



(eie 2 -l)6(A) eia(A)+6(A) 



.ei[eia(A) + 6(A)] 



e 2 b{X) 



6(A), 



(117) 



where the variables a(A) and 6(A) fulfill the affine version of the genus five octic plane curve f|57|) . 

C 2 (A) = (eie 2 - l)[a 4 (A) + a 2 (A)6 2 (A) + 6 4 (A)] 2 - (s 1 s 2 - 2)a 4 (A) + (2 - e x e 2 + £ 2 )6 4 (A) 
+ (2 + el- e 1 e 2 + el) [e 2 {a\X) + a 2 (A)6 2 (A) + 6 4 (A)} + o(A)6(A)] o(A)6(A) 



- £ 2 a 2 (A)6"(A) + 2e 2 a(A)6(A) + 1. 



118) 



Similar results for the special branch are now obtained considering the divisors ( I44l46|l56fl03jl 
They can be solved linearly and the expressions for the weights are: 
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The special branch 



KA) = ; TTTTT, d(X) = ^r, (119) 

v ; exp(^if)a 2 (A) + 6 2 (A)' v ; a(A) V ; 

m = « 2 (A) + exp(±if)6 2 (A) 6 2 (A) 

^ l j a(A)[a 2 (A) + exp(±if)6 2 (A)]' 71 j a(A) ' 1 j 



where the variables a (A) and 6(A) satisfy the affinization of the genus five sextic plane curve (11041) . 



C 2 (A) = [a 2 (A) + exp(±i^)6 2 (A)][a 4 (A) + a 2 (A)6 2 (A) + A 4 a(A)6(A)] + 6 2 (A) - a 2 (A) (121) 

An important property of the above Lax operators is that there exists a special value for the 
spectral parameter A in which they become proportional to the permutator. This turns out to be 
the point Ao on both curves (j!18|12ip in which a(Ao) = 1 and 6(Ao) = 0. In fact, considering the 
expressions for the weights 6(A), d(X), /(A) and g(X) at the value A we obtain^, 

L a , fc (A ) = V a , k = J2 ^ ® ef (122) 

Let us now turn our attention to the R-matrix. After some cumbersome simplifications we find 
that the R-matrix has a universal structure for both branches once we write its elements in terms 
of an enlarged set of weights given by a(A), 6(A), 6(A), d(X) and /(A). The form of the R-matrix 
becomes similar to that of Lax operators, namely 

Ra,/3(A, /i) 





a(A,/i)[eSi ) 


<D efi + e 33 ^ C> 


54?] 


+ b(A^)[eS l ®e<» ) +eg>®e<« ) ] 




+ 




5 + e 22 ^ C 


5 4?] 


+ [e^ 2 ® ^21^ + e 21 ^ ® + e 23 ^ <8> e 32 ^ 


+ e 32' ) ® e 23^ 


+ 


d(A,/i)[eS?C 


g) e 3 2 + e 23 C 


9 4?] 


+ exp(±i^)d(A,^)[eS ) ®eg ) +e£ ) ® 


e!?] 


+ 


f(A^)[eJ?C 


5 e 3 3^ + e 33 £> 


>4?l 


+ g(A, / u)[eg ) ®eg ) ] 




+ 


[a(A,/i) + exp(^i-)f(A,/i)][e§ 


5 ® eg>] + [a(A,/i) +exp(±i^)f(A, / ,)][e£ ) ®eS?] 



(123) 



1 Note that for the main branch the direct limit lim 17(A) is indefinite. This indeterminacy is evaluated by 

A->A 

substituting the expression for 6(A) in Eq. (|117[) and afterwards the high weights powers are reduced with the help 
of the curve constrain (|118[) . After carrying on these simplifications wc find indeed that <?(Ao) = 1 
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where the matrix elements a(A,/i), b(A,/i), b(A,/i) d(A,/i), f(A,/i), and g(A,/i) are given by, 

= a(A) + [6(A)a(/x) - a(A)&(/x)]&(//) 
a(/x) 

= [1 - fe(/x)6fc)]q(A)6(A) - [1 - 6(A)6(A)]a(/i)6(/.) ^ 

a(A)a(/i) 



b(A,/i) = 6(A)a(/x) -a(A)6(/i), (126) 

d(A)a(//)[6(A)o(//) - a(A)6(//)] , . 

1 /(A)aM&M+6(A)[l 1 J 

a(/x)[6(A)o(//) - o(A)6(/x)][/(A)o(/x) - &(A)6(//)] 
1 ' W /(A)a(//)&(//)+KA)[l-6()"W] ' 

g(A,/i) = -d(A,/x) d ^ ^ • (129) 

By construction this R-matrix satisfies the unitarity property which now can be stated as, 

R Qj/3 (A, fj)Rp >a (n, A) = a(A, fi)a(fi, A)I 3 g> I 3 , (130) 

being also a regular matrix at the point A , R aj( g(A , A ) = "Pq,^. 

In addition to that, we have verified by means of the algebraic procedure explained in Appendix 
A that the R-matrix satisfies the famous Yang-Baxter equation, 

Ri2(Ai, A2)Ri3(Ai, As)R23(A2, A3) = R23(A2, A3)Ri3(Ai, As)Ri2(Ai, A2) (131) 

which is sufficient condition for the associativity of the Yang-Baxter algebra (I113p . 

We finally remark that given a solution of the Yang-Baxter equation we can easily produce 
other equivalent multiparametric solutions by means of the so-called twist transformations. These 
are natural symmetries underlying the Yang-Baxter algebra and for a more detailed discussion 
see for example references [21]. It turns out that a special type of twist will be helpful to make 
a correspondence between our respective quantum spin-1 chains and that derived in the work by 
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Alcaraz and Bariev [12]. We find that in our case the most general diagonal twist that is compatible 
with integrability has the following form, 

r 2 r 2 r 4 

G{ti,t 2 ,t 3 ) = diag(l,r 1 ,r 2 |r 2 ,r 3 ,^|r 2 2 ,^,^) (132) 

T 2 T\ T\T~2 

where n, r 2 and T3 are free additional parameters. 

It can be checked easily that the following transformed Lax operator and R-matrix 

0fc,o( r i> 7 "2,T3)L a)Jfc (A)^(ri,r2,r 3 ) and Gp >a {ri, r 2) r 3 )R tt)j9 (A, f^G^in, r 2) r 3 ) (133) 
is still another solution for the Yang-Baxter algebra fU13l) . 

6.2 Spin Chain Hamiltonians 

The expansion of the logarithm of the transfer matrix T(A) around the regular point Ao is 
known to produce a set of mutually commuting operators. Of particular interest is the Hamiltonian 
describing the interaction of nearest neighbors spins variables on the lattice, 

8 N 
H=— lnT(A)| A=Ao = ^H,, fc+1 (134) 

k=l 

where H fc fc+1 = Vk,k+i-§x^k,k+i{^)\)^Xfi and boundary periodic conditions Hjv,iv+i = Hjv,i are as- 
sumed. 

In order to compute the Hamiltonian we just need to take the derivative on the Lax operators 
at the point A = Ao and impose that et(Ao) = 1 and 6(Ao) = 0. The derivatives of the weights a(A) 
and 6(A) are then related under derivation of the algebraic curves (jll8)12ip and evaluating the 
results at A = A . In what follows we shall list the final expressions for the two-body Hamiltonians 
in the Weyl basis: 

• The main branch 
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The spin chain for the main branch is given by, 

tt(±) / c \ t r (fc) (k+1) , (fc) (AH-l)i _ r (fc) 0+1) i (fc) _ £ 2 r (fc) (fc+1) (fc) (fe+l), 

n k,k+l\ £l ' £ 2j — J l[ e ll e ll + e 33 e 33 J L e 21 e 12 ~+~ e 23 e 32 J ~ [ e 12 e 21 ~+~ e 32 e 23 J 

£1 

+ [^f^ + Jxl^e^ + [^ME + h]e^ + + 

+ J 2 exp(±i^) [eg?eg* 1} + e^efA + J 2 exp(^) [eg^ + eg'eg^] 

- [- + h]e^eg +1) (135) 

where the dependence of the couplings Ji and J 2 on the free parameters E\ and e 2 are, 



h = £l ~ £l£2 + 4 and J 2 = V£l£2 ~ 1 . (136) 

We now remark that the one-parameter integrable spin-1 chain found previously by Alcaraz and 
Bariev [12] is in fact a particular case of our Hamiltonian fl 1 3 5 [) when the respective parameters 
sit on the subspace E\ = ±e 2 . Let us denote the two-body Hamiltonian derived in the work [12] 
by Hk,k+i{t P ,t) where e = ± and t p is the single free parameter in the notation of this reference 
0. With help of a special case of the twisted transformation (11321) we have been able to verify the 
following correspondence among two-body operators, 

- H ktk+l (tp, e) = ^=0 fc+ljfc (l, e)H^+i (~7~> ^k+i.ki 1 * e ) + ^ ^ i S k + S k+i] 

where 5| = — eg denotes the azimuthal component of the spin-1 operator. The second term in 
(1137P is proportional to the azimuthal magnetic field and can always be added since it commutes 
with Hamiltonian while the last two terms vanish under periodic boundary condition and do not 
contribute to the volume Hamiltonian. We note from Eq. (1137j) that the twist was only necessary 
to fit the case e = — 1. 

Of course one can use the more general twist transformation (11321) to generate a family of 
exactly solvable multiparametric Hamiltonians. Because this is a diagonal twist it does not spoil 
2 F.C. Alcaraz informed us that the coupling u in the work [T2] should be read as « = y + ^~ £ - > . 
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the U(l) symmetry and the diagonalization of the respective vertex model transfer matrix could in 
principle be tackled by general algebraic framework proposed in |10j . 
• The Special branch 

In this case we have one-parameter spin chain Hamiltonian and the expression for the corre- 
sponding two-body operator is, 

tt(±) /a \ A 4 (fc+1) (fc) (fc+1) , (fc) (fc+1) , (fc) (fc+1) _ (fc) (fc+l)i _ r (fc) (fc+1) , (fc) (fc+1)] 

n fc,fc+ll A 4j ~ a L e ll e ll "r e 33 e 33 + e ll e 33 + e 33 e ll e 22 e 22 I [ e 23 e 21 + e 32 e 12 J 

r (fc) (fc+1) (fc) (fc+l)n _ ( fc ) (fc+1) i_ (fc) (fc+1) , (fc) (fc+1) , (fc) (fc+l)n 

[ e 21 e 12 ~+~ e 23 e 32 J ex Pl =tl g J L e 12 e 32 ~+~ e 21 e 23 + e 12 e 21 + e 32 e 23 J 

(138) 



7 Conclusions 

In this paper we have investigated the Yang-Baxter algebra for three-state vertex model whose 
statistical configurations are invariant by the U(l) invariance but break in an explicit way the parity- 
time reversal symmetry. We argued that the assumption of unitarity of the respective R-matrix 
imposes us that the functional equations derived from the Yang-Baxter algebra are anti-symmetrical 
on the exchange of the Boltzmann weights of distinct Lax operators. This property provides us 
the means to disentangle involved high degree functional relations in a rather systematic way. The 
integrable manifolds are found by intersecting a number of prime divisors associated to polynomial 
equations which are naturally separable on the distinct weights labels. We have been able to 
uncover two families of integrable nineteen vertex models whose weights are lying on bielliptic 
algebraic curves of genus five. For the family having two free parameters this comes about after 
dealing with the problem of the intersection of two projective surfaces: one of them a rational cubic 
surface and the other a cone generated by an elliptic curve. We have pointed out that genus five 
bielliptic curves can generate to standard elliptic curves when the respective free parameters are 
restricted to particular subspaces. 

The Lax operators have a regular point in which they become proportional to the permutator 
and the respective two families of exactly solvable quantum spin-1 chains have been computed. We 
have found that our two-parameter Hamiltonian family generalizes the integrable one-parameter 
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spin-1 chain discovered by Alcaraz and Bariev [12]. We exhibit a relationship between these Hamil- 
tonians when our free parameters are restricted to the subspace £\ = ±£2- We have found that 
the R-matrix has the same general form for both family of vertex models provided we write it as 
function of a suitable subset of Boltzmann weights. The R-matrix is non-additive with respect to 
the spectral parameters and we have verified that it satisfies the Yang-Baxter equation by means 
of computer algebra system. 

A natural question to be asked is whether these integrable vertex models admit an adequate 
description in the framework of quantum groups such as turned out to be the case of chiral Potts 
model [25]. We think that a possible hint into this direction comes from the degeneration of 
the octic plane curve (157|) into a rational curve for specific values of the parameters E\ and e<i 
discussed in Appendix D. This fact has motivated us to search for a relation among the main branch 
Hamiltonian at such particular parameter values and known rational quantum spin-1 chains. To 
this end we have been able to relate the operator H^'^i = 2,^2 = 2) defined by Eq. ( 11351) to 
that of the trigonometric spin-1 chain based on the quantum superalgebra U 5 [Osp(l|2)], often 
referred as the Izergin-Korepin model [26], when the deformation parameter is q = exp(=pi|). 
This suggests that the vertex models obtained in this paper may be also originated from some 
non-generic three-dimensional representation of the U^[Osp(l|2)] superalgebra probably at roots 
of unity. The immediate difficult would be to find the appropriate representation that is able to 
reproduce the pertinent complete intersection the three quadrics by means of the quantum group 
machinery. Hopefully, this observation will prompt further investigations on other mathematical 
properties that are hidden in these vertex models. 

Finally, we expect that the approach used in this paper to solve a number of entangled functional 
equations could also be applied to study the Yang-Baxter algebra associated to generalized nineteen 
vertex models or even high-state vertex models. We plan to investigate some of these problems in 
future works. 
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Appendix A: Extra Functional Relations 

In this appendix we present the remaining functional relations coming from the Yang-Baxter 
algebra (EJ) not presented in the main text. Besides the relations having four terms explicitly 
exhibited in section 14.31 we have the following extra twelve equations, 



sJi a — dc d — f h f 


— hah 


= 0, 


(A.l 


hf a — dc d — haf 


- fh'h" 


= 0, 


(A.2 


hf a — dc d — haf 


- fh'h" 


= o, 


(A.3 


dc d + hh f +fah 


- a/r 


= 0, 


(A.4 


dc d + hh f + fa h 


- a/'/i" 


= 0, 


(A.5 


dc d + hh h — hh h 


— dc d 


= 0, 


(A.6 


ad c — gc d — dh f 


— da h 


= 0, 


(A.7 


ad c — gc d — dh f 


— dah 


= 0, 


(A.8 


cd a —had — fh d 


-dcg 


= 0, 


(A.9 


cd a — ha d — fh d 


-dcg 


= 0, 


(A.10 


cf c — dg d —hcf 


- fc'h" 


= 0, 


(A.ll 


cf c — dg d —hcf 


- ic'h" 


= 0. 


(A.12 



In addition to that we have also functional relations involving five terms. Contrary to what 
happen so far their total number remains unchanged after the solution of the two terms functional 
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relations. Following Table fl2]) we have twelve functional equations which are given by, 



A.13) 
A.14) 
A.15) 
A.16) 
A. 17) 
A.18) 
A.19) 
A.20) 
A.21) 
A.22) 
A.23) 
A.24) 



be b + eg c — dh d — da d — gc g = 0, 

——/—// — / // / // — / „// -/ // 
do b + ha c — gc d — da / — dh h = 0, 

hh d + ia d — bo a + dc (7 — cd h = 0, 

do + ha c — gc a — da / — da a =0, 

be 6 + ca c — dg — fc / — he a =0, 

bo a — gg d — dc f + cd g — dch =0, 

— I II y // / // — / -r/ — 1 11 

dbb + ga c -tea — he a — dg g =0, 
rod — boa + ha a + dc g — ca a = (J, 
be a + eg c — da a — dh d — gc g =0, 

i~ 7' 7" '7" 1 ' r" 7' " i ' 7 " 

bo a — gg d — dc j + cd g — dc h = (J, 

ga c — he a — fc a — do + dg g =0, 

— in 1 u 1 a 1 n 1 u 
be + ca c — dg a — f c / — he a =0. 

An effective way to check that all the above equations are indeed satisfied for the main branch 
is to proceed as follows. After extracting linearly the weights d, f, g, h and h from the divisors 
(I31ll38ll56|l69l) the functional equations become polynomials only in the variables a, 6, b and c for 
both indices labels. In addition to that these weights are constrained by the remaining divisors 
( ]4Tf43l) . The main idea of our procedure is to replace in a given functional equation powers of a 
subset of variables with the help of the last two divisors: 

• Step One 

We have already mentioned that the functional equations depend only on even powers of the 
weight c. The power c 2 can easily be extracted from the divisor (J4TJ) or equivalently from the 

surface flHUl) . Denoting this amplitude by auxc we obtain, 

ab 2 + Eia?b + e 2 b 2 b + ab 2 ,. .. 

auxc = (A. 25) 

e 2 b 

Now we inspect the highest power in a given polynomial equation denoted here generically 
by eq[*]. Assuming that this power is for example six the dependence on the weight c can be 
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systematically replaced using the following Mathematica code, 

eql = Factor[eq[*]] (A.26) 

eq2 = Factor[eql /. {[c'f ->• auxc' [c'] 4 , [c"] 6 -> auxc" [c"] 4 }] (A.27) 

eq3 = Factor[eq2 /. {[c'] 4 ->■ auxc' [c'f, [c"] 4 -> auxc" [c"] 2 }] (A.28) 

eq4 = Factor[eq3 /. {[c ] 2 — )■ auxc , [c ] 2 — >• auxc }], (A. 29) 

where auxc' and auxc" are given by Eq. 11A.25l) with weights labeled by ' and ", respectively. 
• Step Two 

We next use the same method to eliminate other underisable powers now with the help of the 
divisor ( 1431) . For example, we can use this divisor to eliminate the terms that contain powers higher 
or equal to four in the weight b. Denoting the quartic power on b by auxb we find that its expression 
from surface (18T1) is, 

auxb = e\b 2 [e 2 ab + a 2 - e x ah - b 2 } + e x e 2 b 2 [b 2 + &(e x a + 6)] - [b 2 + etab] [b 2 + b{e l a + 26)] 

(A.30) 

Considering that highest power in the resulting polynomial eq4 on the weight b is seven the 
underisable terms can be replaced as follows, 



eq5 = Factor[eq4/. {[&'] 7 ^ auxb [b'} 3 , [6"] 7 ^ auxb ] [b"] 3 }} (A.31) 

eq6 = Factor[eq5 /. {[b'f auxb' [b'] 2 , [b"f auxb" f] 2 }] (A.32) 

eq7 = Factor[eq6 /. {[b'} 5 ->■ auxb' [b'\, [b"} 5 -> auxb" [&"]}] (A.33) 

eq8 = Factor [eq7 /. {[b'} 4 -> auxb', [b"] 4 -> auxb"}], (A. 34) 

where auxb and auxb are obtained from ( 1A.30I) by using the respective label on the weights. 
• Step Three 

It turns out that the final polynomial relation eq8 is either automatically zero or becomes 
proportional to the factor Aq — A + 1. In the latter case we can use the simple substitution, 

eqend = Factor[eq8 /. A 2 A - 1] (A.35) 
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We finally remark that similar verification for the special branch is much simpler since the 
weight 6 can be easily extracted from the divisor ( H9|) . In this case the functional relations become 
dependent only on the weights a, 6 and c. We now use the algebraic plane curve (11041) to extract 
the power c 4 and denoting it by buxc we obtain, 

a 6 + a 4 6 2 + A a 4 6 2 + a 2 b 4 + A a 2 6 4 + A 6 6 + A 4 a 3 6c 2 + A A 4 a6 3 c 2 , 4 , 

buxc = - 9 — A.36 

cr — b z 

We find that the highest power on the weight c is always governed by c 4 and thus a given 
polynomial eq[*] can be verified through the steps, 

eql = Factor[eq[*]] (A.37) 
eq2 = Factor [eql /. {[c'] 4 -»■ buxc', [c"] 4 buxc"}] (A.38) 
eqend = Factor[eq2 /. A 2 -»■ A - 1] (A.39) 

where buxc and buxc are determined in terms of Eq. ( 1A.36I) . 



Appendix B: Elliptic Curves 

In this appendix we shall show that the cone (18 ip is in fact defined over an elliptic curve. To 
this end it is sufficient to work in a given affine chart and here we choose the one defined by setting 
6=1. In this affine chart we can re-scale the coordinates as follows, 

a = xb and 6 = yb, (B.l) 

and by substituting this re-scaling of coordinates in Eq. (I8ip we find that the polynomial ^(xS, yb, 6) C 
C[x, y] becomes, 

S 2 (x, y) = (el - e\y 2 )x 2 + [2 £l - (ele 2 + e\- 2e x - exe\)y 2 \ x 

+ l + (2-e 1 e 2 + e 2 2 )y 2 + (l-e 1 6 2 )y 4 . (B.2) 

We now can complete the square on the variable x making it possible the elimination of the 
linear term on x of the polynomial ( IB. 21) . This is done at the expense of adding an extra factor 
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that depends only on the variable y and together with the last term of Eq. (1B.2l) results in a 
polynomial that factorizes into two pieces. As a result we are able to define the following one-to- 
one transformation, 



~ . , Ri (yi)xi-R2(yi) 

x f Ra(vi) 

y 1 — > yi 

where the expressions of the polynomials Ri(yi), Rzivi) an d Rziyi) are, 



(B-3) 



Ri(yi) = 4ie 2 Vl - s x e 2 {e\ - elyl)yx, 

R 2 ( yi ) = 2i(e?-eiy?)[2e 1 -(e?e 2 + e3-2e 1 -e 1 el)y 1 2 ], (B.4) 



R 3 ( yi ) = ^{ei-elylf 



The map flB.3fB.4j) is known in the literature as de Jonquieres transformation and the zero set 



of S2(x,y) in the new variables x\ and y\ turns out to be equivalent to the curve, 
2 _ 4 [8 + - 2efe 2 + 4ej + e| - 2 gl e 2 (4 + gj) + ef (4 + 3gp] a 

4^1 — £i£ 2 J 

The above polynomial has already the form of an elliptic curve since the right-hand side of 
Eq. flB.5l) is a degree four polynomial in C[?/i]. In fact, this curve can easily be brought into the 
form of a Jacobi quartic. Let ±7*1 and ±r 2 be the roots of the biquadratic polynomial on the variable 
y\. Then by means of the straightforward re-scaling of coordinates, x\ = rir 2 x 2 and y\ = r\y 2 
the plane curve flB.5|) can be rewritten as 



*2 = (1-2/2X1-^2), (B.6) 

whose corresponding modulus parameter is k — ^. 

We remark that all the above reasoning is valid as long as the discriminant of corresponding 
biquadratic polynomial on the variable y\ is not zero. By direct inspection of the right-hand of 
Eq.f ]B.5j) one finds that the expression of such discriminant is, 

A = [8 + e\ - 2e\e 2 + Ae\ + e A 2 - 2e x e 2 (A + e\) + e\{A + 3^)] 2 - 64(1 - e^) 2 , (B.7) 

It turns out that when A = the elliptic plane curve (IB.5P can be factorized in terms of 
two conies and therefore the original surface (jSTj) becomes rational ruled. From Eq.f lB.7p it is not 
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difficult to find this generation occurs in the following one-dimensional submanifolds, 

4 - 2ei + e\ - 2e 2 - e x e 2 + e\ = 0, (B.8) 
A + 2e 1 +e\ + 2e 2 - e x e 2 + e\ = 0, (B.9) 
4e? + e \ - 2e\e 2 + Ae\ + ?>e\e\ - 2e x e\ + 4 = 0. (B.10) 

Note that the above constraints are symmetrical under the exchange of parameters E\ -H- e 2 
emphasizing that our initial choice of free parameters was indeed appropriate. In addition, we 
observe that the first two submanifolds (1B.8HB.9I) can be further reduced as the product of linear 
terms given by Eq.f llOOflOip . respectively. 



Appendix C: Singularities of Curves 

The purpose of this section is to present the technical details concerning the singular locus of 
the degree eight algebraic curves discussed in subsection 15.11 
• The Curve Ci(a, b, c) 

We start by recalling that the singularities of this curve are sited in the following points, 

7T 2lV 

[a s :b s :l], [±exp(i-) : 1 : 0], [±exp(i— ) : 1 : 0], (C.l) 

where the coordinates a s and b s are a subset of solutions of the relations, 

e 2 2 al + (ei£ 2 - 2)a s b 2 s + e 2 b s = 0, (C.2) 
e\bl + (2 - e x e 2 + e 2 2 )a% + e 2 a s = 0. (C.3) 

The above equations can be solved by first considering the resultant of the polynomials with 
respect either to a s or b s . The resultant has the merit to eliminate one of the variables and as result 
we have an univariate polynomial. Considering the resultant with respect to the coordinate a s we 
find, 

b s K + [«2 - eH±a x + et)] bj - a\bl) = 0, (C.4) 
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where the dependence of the coefficients a% and a 2 on the parameters e\ and e 2 are, 

a x = 4 + e 2 [2e 2 + e?e 2 + e\ - e x (4 + el)] , (C.5) 
a 2 = -16 + 24 £l e 2 - 4(1 + 3e?)e£ + 2ei(2 + e?)e^ - (4 + el)e\ + 2e x e\ + e!j. (C.6) 

The trivial solution b s = to Eq. (lC.4j) implies also a s = which has to be discarded since the 
point [0:0:1] does not belong to the curve Ci(a, b, 1). This means that the allowed values for the 
coordinate b s are the eight roots of the second factor of the polynomial (IC.4I) . The corresponding 
coordinates for a s can now be obtained by applying similar reasoning we used to intersect the 
surfaces in subsection 15.11 They can be expressed in terms of the variables b s by the following 
expression, 



£2 



{a 3 + [e 2 (ei - eg) - 2]qj6;} 
ai + £ 2 [2 -e 2 (ei - e 2 )] 



(C.7) 



where the coefficient 0:3 is given by, 



a 3 = -32 + 64eie 2 -8(7 + 6e?)el + 4e 1 (21 + 4e?)e^-2(24 + 21e? + e^)e^ 

+ ei(48 + 7e?)e^-2(ll + 6e?)e^ + lleie^-5e^. (C.8) 

• The Curve Q±(x, y, z) 

By solving the polynomial equations associated to the singular locus ( ITS]) of the target curve of 
the double cover map <j) we find nine singular points. Five of them are located on the affine plane 
and they are given by 

P A = [0 : : 1] and [x s : y s : 1], (C.9) 
where the coordinates x s and y s satisfy the following decoupled equations, 



a\x\ + a 2 x 2 s - e\ = 0, (CIO) 
e 2 (2 - eie 2 ) - e 2 a 1 x 2 s 
4 - 4eie 2 + e\e\ + e\ 



_ t 2 ^ - ziz 2 ) - t2UlA s (n s 

y<* — ~ 7TZ — ~^2J2~T~J- l u - ii J 



The singularities at the infinity line z = sit on the same places of corresponding singular 
points associated to the domain degree eight curve, that is 

7T 27T 

P oo = [±exp(i-):l:0], [± exp(i— ) : 1 : 0], (C.12) 
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except that now they behave as ordinary double points. 

The only singularity that is not an ordinary double point turns out to be the one sited at the 
origin of the affine plane Pq = [0 : : 1]. The respective index of multiplicity is mp =4 and it has 
an extra neighboring infinitesimal singularity. The desingularization diagram is thus given by, 

Qi ^ Qi ^ QiKM (C- 13 ) 

where the curve Qi carries the infinitely near singularity associated to the point Pq = [0 : : 
1] whose index of multiplicity is also four. By using this information we can easily obtain the 
corresponding genus of the normalization Q 1; see Eq. (198ft . 



Appendix D: Reducible Curves 

The purpose of this Appendix is to present the explicit expressions of the plane curves resulting 
from the intersection of the surfaces (I80|l8ip when the parameters E\ and are restricted to the 
submanifolds (I100H102]) . We shall also see such degeneration gives origin to singular algebraic 
quartic curves with genus one. 

• The linear submanifolds: 

We first notice that it is enough to consider one of the submanifolds (HOOflOip because they are 
trivially related under the transformation e\ — > —e% and 82 —> —£2- In the case of the subman- 
ifold ( HOOp we find, after eliminating the parameter s 2 , that the octic plane curve (15T|) becomes 
factorizable in terms of the following product, 

Ci(a, b, c) = g(a, b, c)g(-a, b, ic) (D.l) 
where the expression of the degree four plane curve g(a,b,c) is, 

g(a,b,c) = (~ ±i^ + £l )a 4 ±(-l=Fi>/3-ei± iV3e 1 )a 3 b + ±i^ T V3 £l )aV 
± (-2 + e x ± iv^i)^ 3 + (i T i^ - £l )b 4 ± (-1 ± i\/3 + ex)a 2 c 2 
+ (-2 + |±i^s0afec 2 ±(-l T iV3-|±i^ £ 0feV + (-i±i^)c 4 (D.2) 
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such that the symbol ± refer to the two possible signs of the submanifold fllOOp . 

The plane curve (ID. 20 has only one singular point sited at the infinity line c = which behaves 
as a tacnode and as a consequence of that it has genus one. This means that there exists room for 
emerging new singularities for particular values of the parameter e± and an extra generation to a 
rational plane curve may be possible. Indeed, we find that this happens at the following values of 
the quadratic sub manifolds ( llOOflOip . 

2i . 
E\ = £ 2 = ±2, and E\ = —£ 2 = ±—j=, (D-3) 

v3 

where now the quartic plane curves factorizes once again in the product of two conies. As a result, 
the original octic plane curve AST]) degenerates to genus zero curve and the corresponding weights 
of the vertex model becomes trigonometric at the specific points (1D.3j) . 
• The quartic submanifold: 

The essential ingredient in the analysis of the third submanifold (11021) is to note that we are 
in fact dealing with a curve that can be rationally parameterized. In fact, one easily finds that 
submanifold (11021) has the maximal number of three ordinary singular points when viewed in the 
projective space. Now given such three double points on a quartic curve we can fix any other point 
on it and pass conies through these four points. By using this family of conies it is well known that 
one can establish a rational parameterization and the final result is, 

4(-3 + 14k) (147t 2 - m + 294/ct + 12 - 46k) (-49t 2 + 98nt + 2 - 24k) 
£l (133t 2 - 56t + 84Kt + 4 - 6/«)(343t 2 - 196* + 980/rf + 24 - 190/?) ^ ' ' 



e-2 



2(-4 + 63K)(147t 2 - 98* + 294k* + 12 - 46K)(637t 2 - 294t + A90nt + 36 - 138k) 
13(133t 2 - 56* + 84Kt + 4 - 6K)(343t 2 - 196* + 980k£ + 24 - 190k) 

(D.5) 



where t is the parameterization variable and k is a constant factor k = 6± ]^ ■ 

With the help of the parameterization (ID.4j) we then are able to investigate the explicit factor- 



ization of the octic plane curve (IBTj) in terms of the product of two quartic curves, namely 

Ci(a, b, c) = h(a, b, c)h(—a, b, ic) (D.6) 



50 



where the expression for h(a, b, c) is, 



Kl , 4 , .4 , 2^ , (29 ± 36iV^3> 2 4 2(±2iVl3 + 9)«1 

(a + o + a b H c + 

109531219 v ; 14567652127 6900466797 



2(+2iVl3 + 9)^ 2 2 2(±9i7l3-26)K 3 / t4 2 
1684480617001 ° 388726296231 ° ° 1 ' j 

such that the coefficients «i, • • • , K4 are given by, 
k x = (45619t 2 - 16856i ± 392i v / l3t + 1836 + 86iv / 13)(160 - 1372* + 2401t 2 + 6i>/l3) 

*2 1 non:. /To-/- I7n^ -■- i nn:. /To , o^hcn on^m±o/:.& , cciw2 



k 2 = (16807r ± 980iV13t - 3724t + 190iV13 + 36) (160 - 2240t ± 84iV13t + 6517r + 6iV13) 
k 3 = (147t±3iv / T3-31 + 12i + 2v / 13)(147t-31±12i-2v / 13±3iv / 13) 
k 4 = (637t±5iv / T3-117 + 26i-12v / 13)(637t±5iv / T3-117±26i + 12\/i3) (D.8) 

For generic values of the free variable t one finds that the quartic plane curve ( ID. 71) has two 
ordinary singular points and thus has again genus one. We remark that in this submanifold we 
have not been able to find a further generation to rational curves. This however can not be ruled 
out since in this case the analysis is more subtle. 
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